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The comment by M. Pitkin [Tl on our EPL article 
“Measurements of Newton’s gravitational constant and 
the length of day” claims to provide evidence that a con¬ 
stant G measurement model with an additional Gaussian 
noise term is “hugely favoured” over models employing 
sinusoidal terms when using a Bayesian model selection 
procedure. Unfortunately, we were unable to replicate his 
claims with our own independent analysis testing the hy¬ 
potheses of the following three scenarios for the G mea¬ 
surements: 

1. Constant value: G = a$. 

2. Constant plus a sinusoidal term with period of ap¬ 
proximately 6 years: G = ao + ai cos + 

61 sin j . 



Fig. 1: Probability density function of the normalized residuals 
about a weighted mean G model (hypothesis 1). 


normalized a values of 4.0, 2.6 and 2.0 for the weighted 
residuals of scenarios 1, 2 and 3, respectively, suggesting 
the two-period model is favored. We also computed his¬ 
tograms for the 19 weighted residuals and best fit proba¬ 
bility density functions for the three hypotheses. Impor¬ 
tantly, the probability density function of the weighted 
residuals about a mean value of the G measurements (hy¬ 
pothesis 1) follows more of a uniform distribution whereas 
for the two-period sinusoidal model (hypothesis 3) the 
probability density function of the weighted residuals ap¬ 
pears to follow a normal distribution, suggesting a pos¬ 
sible error in Pitkin’s analysis. See Figs. [I][3] below for 
our outputs from Mathematica. Thus, we stand by our 


3. Constant plus sinusoidal terms with two different pe¬ 
riods of approximately 6 years and 1 year: G = 
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We used a non-linear regression analysis with a minimiza¬ 
tion of the LI norm to determine the best fit values for 
the input parameters of each of the above cases (i.e. for 
ai, bi and Pi values). After fitting to the G data, we found 
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Fig. 2: Probability density function of the normalized residuals 
about a constant plus two-period sinusoid model (hypothesis 
3). 


conclusions of potential periodic terms in the reported G 
measurements (see our added appendix of [2] in response 
to [3] for our logic with a two-period sinusoidal model). 
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Normalized Residual 


Fig. 3: Probability density functions of the normalized residu¬ 
als for all three hypotheses (red: hypothesis 1, blue: hypothesis 
2, green: hypothesis 3). Overlayed are fitted normal distribu¬ 
tions for each hypothesis with clear indication that hypothesis 
3 residuals have a tighter Gaussian fit than the others. 
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